home *** CD-ROM | disk | FTP | other *** search
/ HyperLib 1997 Winter - Disc 1 / HYPERLIB-1997-Winter-CD1.ISO.7z / HYPERLIB-1997-Winter-CD1.ISO / オンラインウェア / MLT / XTRACE311.sit / XTRACE311 / X-TRACE 3.1.1.rsrc / TEXT_1002_psph.xti.txt < prev    next >
Text File  |  1996-09-01  |  833b  |  61 lines

  1. xtrace    "poly-sphere"
  2. {
  3.     procedure    psph    l,m,n, r, pn,pn2
  4.     {
  5.         nvar    pn=16;
  6.         nvar    pn2=int (pn/2);
  7.         nvar    h=r;
  8.  
  9.         var        t1,t2;
  10.         var        c1,c2,s1,s2;
  11.         var        r1,r2,z1,z2;
  12.         var        x1,y1,x2,y2,x3,y3,x4,y4;
  13.  
  14.         para    l,m,n;
  15.         sim        r,r,h;
  16.         boundary    sphere    0,0,0, 1
  17.         {
  18.             var        i,j;
  19.             j=0;
  20.             while    lt (j - pn2)
  21.             {
  22.                 t1=PI*j/pn2;
  23.                 t2=PI*(j + 1)/pn2;
  24.                 z1=cos (t1);
  25.                 r1=sin (t1);
  26.                 z2=cos (t2);
  27.                 r2=sin (t2);
  28.  
  29.                 i=0;
  30.                 while    lt (i - pn)
  31.                 {
  32.                     t1=2*PI*i/pn;
  33.                     t2=2*PI*(i + 1)/pn;
  34.                     c1=cos (t1);
  35.                     s1=sin (t1);
  36.                     c2=cos (t2);
  37.                     s2=sin (t2);
  38.                     x1=r1*c1;
  39.                     y1=r1*s1;
  40.                     x2=r1*c2;
  41.                     y2=r1*s2;
  42.                     x3=r2*c1;
  43.                     y3=r2*s1;
  44.                     x4=r2*c2;
  45.                     y4=r2*s2;
  46.                     poly
  47.                     {
  48.                         pt    x1,y1,z1;
  49.                         pt    x2,y2,z1;
  50.                         pt    x4,y4,z2;
  51.                         pt    x3,y3,z2;
  52.                     }
  53.                     i=i + 1;
  54.                 }
  55.                 j=j + 1;
  56.             }
  57.         }
  58.     }
  59. }
  60.  
  61.